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Abstract — However evidence was collected which shows that, 
with appropriate design of walls materials can contribute to 
reducing loads and achieving energy savings in buildings, while 
securing a comfortable indoor environment. Among the key 
factors of this design, we have developed a new equation for 
Controller the average Nusselt number and the exchange 
coefficient h to be obtained before practice and actual 
realization of the walls which gives us a prediction before 
conception walls. First, CFD software is used as a technique to 
modelling the behaviour of fluid and the thermal convection in 
the house with different Rayleigh numbers, in this context, an 
analogy was used for the functions which are discretized by the 
finite difference method and integrated in the CFD code which 
is based on the finite volume method. Secondly, we vary the 
thickness of the building material of the outer wall four times 
and calculate the Nusselt number and exchange coefficient of 
heat transfer aims to find a cloud point respectively for the 
thicknesses e = 0, L / 40, L / 20 and L / 10. After, we developed a 
relationship that helps us to know the exchange ratio for each 
thickness (e) belongs to the interval [0, L / 10] by the Lagrange 
polynomial for Rayleigh number equal 10 3 . This method for 
predicting exchange coefficient of convection for to optimize the 
design of walls in buildings. Finally we developed FORTRAN 
program under code blocks basis to facilitate the task of the 
control of this nonlinear equation, in all the point 
[0< e < L/10]. 

Index Terms — Thermal comfort, Building material, Design of 
walls, Natural convection, Polynomial interpolation method. 
Heat transfer. 


I. INTRODUCTION 

Actually, In Algeria, The construction industry must start 
to take into account the changing needs of the population in 
its projects. It’s not just homes that need to be adapted to 
ensure a safe and inclusive environment, but also the wider 
surroundings. Local shops and town centres have to be 
accessible and to accommodate the needs of people. 


Elhadj. Benachour, Departement Of Sciences, Energarid Laboratory, 
Faculty of Science , University of Tahri Mohamed, Bechar, 08000, Algeria 

Belkcem . Draoui, Department of Technology, Energarid Laboratory, 
Faculty of Technology , University of Tahri Mohamed, Bechar, 08000, 
Algeria. ( e-mail: bdraoui@yahoo.com ) 

Bachir . Imine , LASP Laboratory, Faculty of Mechanical Engineering, 
USTO,31000 

Khadidja Asnoune , Department of Technology, Energarid Laboratory, 
Faculty of Technology , University of Tahri Mohamed, Bechar, 08000, 
Algeria 

Mohamed Hasnat , Department of Technology, Energarid Laboratory, 
Faculty of Technology , University of Tahri Mohamed, Bechar, 08000, 
Algeria 

Rahmani. Lakhdar, Department of Technology, Energarid Laboratory, 
Faculty of Technology , University of Tahri Mohamed, Bechar, 08000, 
Algeria 


This survey indicates that the biggest improvements could be 
made to the building sector, which must be properly adapted 
to ensure a safe environment for a workforce .The building 
sector today is known to be consuming 40% of the world 
energy [1,2], and in turn, supports 23-40% of the world's 
greenhouse gas emission, particularly C02 [3]. It has become 
common in recent years to gather data on human attitudes and 
behaviour in building energy research through interviews 
with building occupants and other relevant actors. Examples 
of this kind of research are found in the main journals which 
deal with technical aspects of energy in buildings, including 
Building and Environment [4—12], Building and 
Environment [13-15], Energy Storage [16-20], Energy and 
Buildings [21-33]. In most situations, the mechanical cooling 
devices offer solutions that are neither environment friendly 
nor energy sustainable. The mechanical devices are 
non-functional and cannot offer thermal comfort without 
energy input. Hence utilization of advanced building 
materials and passive technologies in buildings may offer the 
solution for thermal comfort demands, substantially reduce 
the energy demand, impact on the environment and carbon 
footprint of building stock worldwide [34] .Thermal comfort 
is dependent and influenced by a range of environmental 
factors viz. air temperature, radiant temperature, humidity, 
air movement, metabolic rate or human activity [35,36]. 
Also, numerous studies across the world have shown the 
impacts of hot working environments on the working 
population [37—43]. The choice of which material or 
combination of materials is used depends on a wide variety of 
factors. Some insulation materials have health risks, some so 
significant the materials are no longer allowed to be used but 
remain in use in some older buildings such as asbestos fibers 
and urea.Also, the cost can be high compared to the 
traditional insulation. The interior or external insulation does 
not often becomes the only possible solution in particular 
with the old buildings in the Sahara of Algeria. Hygrothermal 
parameters of the existing wall material should be reported 
and well known for designing powerful and durable walls in 
the time study. First, in this paper, CFD software is used as a 
technique to modelling the behaviour of fluid and the thermal 
convection in the external wall of the house with different 
Rayleigh numbers [l.E3< Ra<l.E5], In this context, an 
analogy was used for the functions which are discretized by 
the finite difference method and integrated in the CFD code 
which is based on the finite volume method. Secondly, The 
most important part in this work is to vary the thickness of the 
building material of the outer wall four times and calculate 
the Nusselt number and exchange coefficient of heat transfer 
aims to find a cloud point respectively for the thicknesses e = 
0, L / 40, L / 20 and L/10. After, we developed a relationship 
that helps us to know the exchange ratio for each thickness e 
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belongs to the interval [0, L / 10] by the Lagrange polynomial 
interpolation method, and then we developed a FORTRAN 
program to control the nonlinear equation of order three. This 
method for predicting exchange coefficient of convection to 
optimize the design of walls in buildings before starting the 
wall construction for Rayleigh number equal 10 3 . 


Y-momentum. 


e_v +u w +v d J- 

8t 8X 8Y 


+ ttVV+g/? — 

pdx 8X 
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II. Geometric configuration Energy. 




Adiabatic wall 


Fig. 1. Schematic of the studied configuration 
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The derived equation of motion Eq. (2) over Y and the 
equation of motion Eq. (3) by contributing to x, then, after 
subtracting the two equations obtained, we obtain the 
equations dimensionless variables in writing Helmotz in 
terms of vorticity and stream function formulation are as 
follows: 


^ + U — + V — = P,V 2 l 
dt 8X dY 


f Ra Pr — (5) 

dX 


To simplify the problem, assume that: 

• The fluid is Newtonian and incompressible. 

• The heat dissipation by viscous friction is neglected. 

• The Boussinesq approximation is considered 
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III. Objectives 

We can express our aim for this study in the following 
points: 

• Prediction of the wall design through numerical 
evaluation of the convection in buildings by the 
Lagrange polynomial interpolation method. 

• Modelling Study of the convection - conduction 
coupling. 

• Study of the effect of the distribution of the heat 
inside buildings on the convection. 

• Comparative study between the construction of 
buildings with different thicknesses and without 
insulation. 

• Study the effect of the Rayleigh number on the 
convection in building. 


The stream function and vorticity are related to the velocity 
components by the following expressions: 
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The dimensionless parameters in the equations above are 
defined as follows: 



(9) 


IV. Mathematical model 
The fluid is assumed incompressible and obeys the 
Boussinesq approximation. In these cases, continuity in two 
dimensions and the equations governing the flow and energy 
is given by: 


For solid (Concrete), we are interested only in the following 
heat equation: 
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Continuity. 


dU 

dX 


l 5V = 0 
dY 


( 1 ) 


X-momentum. 


™ + U™ + V™ 
dt dX dY 
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V. Procedure of Simulation 
First, The numerical calculation was conducted using the 
GAMBIT and FLUENT software. The numerical procedure 
used in this work is that of finite volume. It involves the 
integration of differential equations of mathematical model 
on finite control volumes for the corresponding algebraic 
equations. In this context, an analogy was used for the 
functions which are discretized by the finite difference 
method and integrated in the CFD code which is based on the 
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finite volume method. The SIMPLE algorithm [44] was 
chosen for the coupling speed pressure in the Navier-Stokes 
equations on a staggered grid. The convective terms in all 
equations are evaluated using the schema apwind first order. 
The discretization of the time term is made in a totally 
implicit scheme. The convergence of the solution is 
considered reached when the maximum relative change of all 
variables (u, v, w, p, t) between two successive time is not 
less than 1 ,E4. With an aim of following well any variation of 
the fields thermal and hydrodynamic, we used a uniform grid 
of 14241 nodes and 14480 Elements in non-stationary mode. 
Secondly, we developed a FORTRAN program to control the 
nonlinear equation of order three. This method of 
interpolation was used for predicting exchange coefficient of 
convection to optimize the design of walls in buildings 

VI. Validation of the Code of Calculation 


to = 0 if i*j 
[W = 1 if i = J 

7.1 The polynomial interpolation points 


We wish to find the polynomial interpolating the points 


Table2 The sets of polynomial interpolation points. 


Nu a 


0 

2.24 


1/40 

1.59 


1/20 1/10 
T35 098" 


Where, the points of Nu (e) are obtained by numerical simulation 
of the FLUENT software, for the Rayleigh number equal 10 3 . 


VI. 1 Validation of the model 

The computer code was validated on a natural convection 
problem of stale air in a square cavity with vertical walls 
differentially heated and adiabatic horizontal walls. Our 
results were compared with those obtained by De Vahl Davis 
(1983) [45] .The latter dealt with the same problem by 
adopting the finite difference method with the 
vorticity- stream function formulation (see Table 1). 


Table 1 Comparisons of the results of validation 



Ra = 10 3 

Ra = 10 4 

Ra = 10 3 

Nu 

present 

present 

Present 

(average 

study: 

study : 

study : 

number) 




(hot wall) 

1.2 

2.257 

4.64 

Nu 

De Vahl 

De Vahl 

De Vahl 

(average 

Davis 

Davis 

Davis 

number) 

[45] 

[45]: 

[45]: 

(hot wall) 





1.12 

2.243 

4.52 


VII. Mathematical model for The Lagrange 

POLYNOMIAL INTERPOLATION METHOD 
The Lagrange interpolating polynomial is 
the polynomial of degree <(n-l) that passes through the n 
points (xl, yl=f(xl)) , (x2, y2=f(x2)) , ..., (xn, yn=f(xn)), 
the interpolation polynomial in the Lagrange form is a linear 
combination and is given by : 

p n = W (11) 

For Lagrange basis polynomials 

fl(x-Xj) 

)= o 

L i ix)= J ^ 

no*-*,) 

/=o 

Where 


VIII. Results and Discussion 

The boundary conditions have been established to simulate a 
geometric configuration used frequently in two-dimensional 
approximation .The structure of the flow, the temperature 
field and heat transfer through the hot wall are discussed in 
this section.In this study, to target the most important goal, 
we will show and studied the dynamic and thermal behavior 
of the fluid in the cavity. 

VIII. 1 Isotherms 

The isotherms are shown in Figures [2], The heat 
distribution in the cavity is in accordance with the fluid 
circulation revealed by isotherms and iso-currents. Indeed we 
find a heating fluid from the interface, if it causes the change 
of the heat distribution in the cavity (see figure) for different 
numbers of Rayleigh For a fixed value of the number of 
Prandtl equal 0.71 (air fluid). Gradually, as the Rayleigh 
number has increased, the isotherms become increasingly 
wavy and heat transfer increases, so the flow intensifies and 
natural convection is expanding and predominates (natural 
convection is predominant). 



Fig. 2. The Isotherms for wall thickness e = L/20, and for 
different Rayleigh Numbers, Pr=0.71. 


n 


Kx-Xj) 
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VIII. 2 Lagrange polynomials 

The points of Nu (e) are obtained by numerical simulation of 
the FLUENT software, for the Rayleigh number equal 10 3 . 
(See table 3) 


1120 3 68 2 7.066 , 

PJx) = -x 3 + — X 2 x+1.12 

3 3 L 3 L 2 L 

Therefore, we can write, 

— , 1120 3 68 2 7.066 , .. 

Nu(e) = —e+—e e+1.12 


(18) 


(19) 


Table3 Nusselt number for different values of thickness, and 
for Ra=10 3 


Ra 

1 

Xi=e 

Yi=Nu (e) 
Average Nusselt 
number 

10 4 

0 

0 

1.12 

10 4 

1 

e=L/40 

0.98 

10 4 

2 

e=L/20 

0.89 

10 4 

3 

e=L/10 

0.72 


If (0, 1.12),(L/40, 0.98), (L/20, 0.89), (L/10, 0.72), are given 
data points, then the cubic polynomial passing through these 
points can be expressed as, 


P 3 (x) = y 0 l 0 + y l l 1 +y 2 l 2 + y 3 l 3 


We would have the four basis polynomials 

lo J X -AO>^-L^ X -to> 

(O-^A (O-^-) (O--^) 


Nu : Average Nusselt number, 
e : Wall thickness. 

Also, we can write: 

Nu(e) = -373.333 5 3 + 68 B 2 -7.066 5+1.12 (20) 

Where: B: e/L 
Therefore, we can write. 


In heat transfer within a fluid, the Nusselt number (Nu) is 
the ratio of convective to conductive heat transfer. The 
average value of the heat transfer coefficient is often needed 
in process engineering and design applications. In this 
context, the average Nusselt number is given by: 


Nu 


hH 

A 


=> h = 


ANu 

H 


( 21 ) 


Where h is the convective heat transfer coefficient of the 
flow, H is the characteristic length, X is the thermal conductivity of 
the fluid. 


VIII.4 Average Nusselt numbers 


O-O) O-— ) (x — 


40 20 40 


«*- ° )c ?:.4o )c * .o> 

io- 0) ^-io 5 




We can see that each term is of degree three, so the entire 
sum has degree at most three. 

The polynomial P(x) given by the above formula is called 
Lagrange’s interpolating polynomial and the functions Eq. 
(14, Eq). (15), Eq). (16), Eq.( 17) are called Lagrange’s 
interpolating basis functions. 

VIII. 3 Nusselt number correlations 

It should be noted that the numerical result given by the 
equation Eq.(20). Remarkably similar to the estimate of the 
average Nusselt number by the Lagrange interpolation 
method for each wall thickness between [0, L / 10] but only 
forRa=10 3 . 


The influence of the Rayleigh number and Wall thickness on 
the average Nusselt number is shown in Fig. 5 and fig 6. The 
Nusselt number decreases with increasing of the wall 
thickness. Therefore, the heat transfer also decreases because 
the inertia and the thermal resistance of the wall increases. 
The average value of the Nusselt number increases as the 
value of the Rayleigh number (Ra) increases. Broadly, 
advection becomes stronger and thus heat transfer increases. 
This trend is qualitatively consistent with that seen in the 
convection for important values of Ra. 



Fig. 5. Profiles of average Nusselt number for different wall 
thickness and for different Rayleigh number, Ra=10 3 ,10 4 ,10 5 , Pr 
= 0.71. 
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Fig. 6. Profiles of average Nusselt number for different Rayleigh 
number, Ra=10 3 ,10 4 ,10 s , Pr = 0.71. 

IX. Conclusion 

Extensive numerical results for the laminar natural 
convection heat transfer from a heated wall are obtained for 
the range of conditions as: Rayleigh number, 10 3 < Ra <10 5 
and Prandtl number, Pr =0.71, Concrete wall for different 
thicknesses is viewed with 0< e < L/10. We are interested in 
convection conduction coupling. First, CFD software is used 
as a technique to modelling the behaviour of fluid and the 
thermal convection in the externel walof the house, in this 
context, an analogy was used for the functions which are 
discretized by the finite difference method and integrated in 
the CFD code which is based on the finite volume method. 
Secondly, The most important part in this work is to vary the 
thickness of the building material of the outer wall four times 
and calculate the Nusselt number and exchange coefficient of 
heat transfer aims to find a cloud point respectively for the 
thicknesses e = 0, L / 40, L / 20 and L / 10. After, we 
developed a relationship that helps us to know the Nusselt 
number and exchange ratio for each thickness ( e ) belongs to 
the interval [0, L / 10] by the Lagrange polynomial 
interpolation method, and then we developed a FORTRAN 
program to control the nonlinear equation of order three 
(equation 19. or. 20). This method for predicting exchange 
coefficient of convection for to optimize the design of walls in 
buildings. Since, Interpolation is the process of defining a 
function that takes on specified values at specified points. 


Predictive simple correlation was developed to estimate the 
value of the average Nusselt number and the coefficient of 
heat transfer for exchange any thickness ranging from 0 to 
L/10 .it just enough to replace the thickness value in equation 
19, to calculate the Nusselt number planned before the design 
and construction of walls 
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